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ABSTRACT 

We have exploited the new, deep, near-infrared UltraVISTA imaging of the COSMOS field, 
in tandem with deep optical and mid-infrared imaging, to conduct a new search for luminous 
galaxies at redshifts z ~ 7. The year-one UltraVISTA data provide contiguous Y, J, H, K s 
imaging over 1.5 deg 2 , reaching a 5er detection limit of Y+ J ~ 25 (AB mag, 2-arcsec di- 
ameter aperture). The central ~ 1 deg 2 of this imaging coincides with the final deep optical 
(u*,g,r,i) data provided by the Canada France Hawaii Telescope (CFHT) Legacy Survey 
and new deep Subaru Suprime-Cam z'-band imaging obtained specifically to enable full ex- 
ploitation of UltraVISTA. It also lies within the Hubble Space Telescope (HST) / 81 4-band 
and Spitzer IRAC imaging obtained as part of the COSMOS survey. We have utilised this 
unique multi-wavelength dataset to select galaxy candidates at redshifts z > 6.5 by searching 
first for Y+ J-detected objects which are undetected in the CFHT and HST optical data. This 
sample was then refined using a photometric redshift fitting code, enabling the rejection of 
lower-redshift galaxy contaminants and cool galactic M, L, T dwarf stars. The final result of 
this process is a small sample of (at most) ten credible galaxy candidates at z > 6.5 (from 
over 200,000 galaxies detected in the year-one UltraVISTA data) which we present in this pa- 
per. The first four of these appear to be robust galaxies at z > 6.5, and fitting to their stacked 
spectral energy distribution yields z p hot = 6.98 ± 0.05 with a stellar mass ~ 5 x 10 9 M Q 
and rest-frame UV spectral slope /3 ~ —2.0 ± 0.2 (where fx oc A^). The next three are also 
good candidates for z > 6.5 galaxies, but the possibility that they are dwarf stars cannot be 
completely excluded. Our final subset of three additional candidates is afflicted not only by 
potential dwarf-star contamination, but also contains objects likely to lie at redshifts just be- 
low z = 6.5. We show that the three even-brighter z > 7 galaxy candidates reported in the 



COSMOS field by Cap aket alT] < |201 1] > are in fact all lower-redshift galaxies at z ~ 1.5 - 3.5. 
Consequently the new z ~ 7 galaxies reported here are the first credible z ~ 7 Lyman-break 
galaxies discovered in the COSMOS field and, as the most UV-luminous discovered to date at 
these redshifts, are prime targets for deep follow-up spectroscopy. We explore their physical 
properties, and briefly consider the implications of their inferred number density for the form 
of the galaxy luminosity function at z ~ 7. 

Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift. 



1 INTRODUCTION 

The advent of deep near-infrared imaging on the Hubble Space 
Telescope (HST) ushered in a new era in the discovery of galax- 
* E-mail- raab@roe ac uk i es at redshifts z ~ 7 — 8. Installed in 2009, the infrared channel of 

f Scottish Universities Physics Alliance Wide Field Camera 3 (WFC3/IR) has already been used to obtain 
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near-infrared imaging to detection limits Y, J,H ~ 29 (AB mag). 
This ultra-deep near-infrared imaging, when combined with exist- 
ing deep optical imaging in fields such as the Hubble Ultra Deep 
Field (HUDF) and GOODS-South, has enabled the first detection 
of significant numbers of Lyman-break galaxies (LBGs) at z > 6.5, 
and hence the first meaningful studies of the rest-frame ultraviolet 
(UV; A r est — 1500 A) galaxy luminosity function (LF) at z ~ 7 
(e.g.|McLure et al.|2010||Oesch et al.|2010[|Bouwens et al.|2010l 
|Finkelstein et al.|2010||McLure et al.|201 l||Bouwens et al.|201 1| >. 

The depth of this HST imaging, coupled with the (relatively) 
small field-of-view of WFC3/IR has meant that, to date, these new 
studies of the z ~ 7 LF have been largely focussed on the faint 
end, with the WFC3/IR samples dominated by sub-L* galaxies 
(Misoo — "20 — > —18 mag). This has of course been of enormous 
value, as it has revealed that the faint-end slope of the z ~ 7 LF is 
steep (e.g. a = -1.7, |McLure et al.|2010[ a = -2.0, |Bouwens| 
|et al.|20l"l) , implying that it is the large population of fainter galax- 
ies which is likely responsible for the re-ionisation of the Universe 
(e.g. |Robertson et al.|2010| l. Indeed, to further explore the faint end 
of the z ~ 7 LF, and extend such studies out to z ~ 9, even deeper 
WFC3/IR imaging of the HUDF will be completed before the end 
of 2012 (HST Cycle- 19 program GO 12498). 

However, at the same time the importance of better determin- 
ing the number density of brighter galaxies at this crucial epoch 
has not been over-looked. In particular, the 3-year, 900-orbit Cos- 
mic Assembly Near-infrared Deep Extragalactic Survey (CAN- 
DELS) HST Treasury Program ( Gro gin et al.||20lT] |Koekemoer| 
|et al.|[20TT| >, reaching detection limits of J, H ~ 27 mag, is al- 
ready providing better statistics on the number density of z ~ 7 
galaxies at luminosities around the (apparent) break luminosity of 
the LF (i.e. at Mf 500 ~ -20.1; |Ouchi et al.|2009||McLure et aT] 
[20T0{|Bouwens et al.|20lT] >. Meanwhile, CANDELS and parallel 
HST WFC3/IR imaging programs such as Brightest of the Reion- 
izing Galaxies survey (BoRG, Trenti et al. 20iT| are also now 
yielding the first significant numbers of brighter galaxies at z ~ 8 
(e.g. |Oesch et al.|2012|[Bradley et al.|20T2l >. 

However even CANDELS, when complete, will cover an area 
of only ~ 800 arcmin 2 . There thus remains a key role for wider- 
area, albeit shallower, ground-based near-infrared imaging to bet- 
ter constrain the number density of rarer more luminous galaxies 
at z > 6.5, and hence properly determining the bright end of the 
galaxy LF at these epochs. The value of degree-scale imaging sur- 
veys of the high-redshift Universe (sampling com oving volumes 



McLure et al. 



~ 100 x 100 x lOOMpC 3 ) was demonstrated by_ 
p006l[2009l l and |Curtis-Lake et al.| ( [2Tjl"2"l >, who utilised the Subaru 
Suprime-Cam and UKIRT WFCAM imaging of the SXDS field to 
search for rare bright galaxies at z ~ 5 and z ~ 6, complementing 
deeper smaller-area imaging work with HST ACS (e.g. Bouwcns 
|et al.|20 07 ), and ultimately yielding determinations of the z ~ 5 
and z ~ 6 LF spanning a dynamic range of over five magnitudes. 

An attempt to push ground-based studies of LBGs out to 
z ~ 7 has been made by |Ouchi et al.|(2b09| > and Castellan o et al.| 
(2010a b). However the latter study (with Hawk-I on ESO's Very 
Large Telescope; VLT) has covered an area smaller than the fi- 
nal CANDELS area, while the former was based on pushing the 
silicon-based CCDs in Subaru Suprime-CAM to their red limit, and 
lacked the longer-wavelength (A > 1 /im) near-infrared imaging 
required to confirm a blue spectral slope longward of the putative 
Lyman break (resulting in highly-contaminated galaxy samples). 

Now, however, degree-scale near-infrared imaging reaching 
the depth required to uncover credible galaxy candidates at z > 6.5 
has finally arrived with the Ultra VISTA survey. This 5-year pub- 
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Figure 1. The multi-band coverage map of the UltraVISTA/COSMOS field 
utilised in this study. Working from the outside in, the year-one 1 .5 deg 2 Ul- 
traVISTA imaging is shown as the large red rectangle, with the HST/ ACS 
7814-band coverage indicated by the jagged outline. The blue irregular 
shape within the UltraVISTA data is the Subaru z'-band mosaic, formed 
from four individual Suprime-Cam pointings. Finally, the central green out- 
lined square is the CFHTLS D2 optical data. This final square area, cover- 
ing ~ 1 deg 2 , is the area utilised in this study, as it contains all the required 
overlapping multi-wavelength data. 



lie survey on the COSMOS field (|Scoville et al.||2007b) with the 
near-infrared camera VIRCAM (Dalt on et al.||2006} on the new 
VISTA survey telescope at Paranal in Chile (Emerson & Suther-| 
|land||20f0| > commenced in 2010, and the first public data re- 
lease of the fully-reduced year-1 d ata was made through ESO in 
Feb 20lf]McCr acken et al.|2012 i. This new near-infrared imag- 
ing covers 1.5 deg 2 to depths of Y = 24.7, J = 24.5, H = 
24.0, K s — 23.8 (5a; 2-arcsec aperture diameter). At the longer 
wavelengths this is not yet as deep as the 0.8 deg 2 imaging being 
delivered in the SXDS field by the ongoing UKIDSS Ultra Deep 
Survey (UDS; Lawrence et al.|200"7] >. However, as well as cover- 
ing a larger area, UltraVISTA provides the crucial new ingredient 
of the first large-area Y-band imaging reaching ~ 25 mag (due, 
in part, to the vastly superior short-wavelength sensitivity of the 
VISTA Raytheon detectors). As discussed further below, this deep 
Y-band imaging is of vital importance, both for the effective selec- 
tion of credible LBGs at z > 6.5 and for the robust rejection of cool 
brown-dwarf star contaminants. Also crucial is the new availability 
of very deep Subaru z'-band imaging over the central 1 deg 2 of the 
UltraVISTA field, obtained over the last three years with the refur- 
bished Suprime-Cam (equipped with its new red-sensitive Hama- 
matsu CCDs); these data were obtained specifically to provide red 
optical imaging of the depth necessary to exploit UltraVISTA in the 
search for extreme-redshift LBGs. 

We note that UltraVISTA does not represent the first near- 
infrared imaging of the COSMOS field. Indeed, [Capak et al.|(20TT] l 
exploited the somewhat shallower H, K s imaging obtained with 
WIRCam on the Canada-France-Hawaii Telescope (CFHT) and 



1 |http : //www. eso ■ org/sci/observing/phase3/data_ 
|releases ■ html] 
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Table 1. The 5cr limiting magnitudes for the relevant optical and near- 
infrared data used in this study, obtained from the rms of the flux from 
apertures placed in blank regions of each image (see Section [3~T) . All mag- 
nitudes were calculated within a 2-arcsec diameter aperture, apart from the 
ACS 7814 values, which used a 0.6-arcsec diameter aperture and the IRAC 
3.6 fim and 4.5 /im values, which were calculated in a 2.8-arcsec diameter 
aperture. 
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26.9 


CFHT/MegaCam 
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27.0 


CFHT/MegaCam 
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26.6 


CFHT/MegaCam 
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26.4 


CFHT/MegaCam 


^814 


26.7 


HST/ACS 
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25.2 


CFHT/MegaCam 


z' 


26.3 


Subaru/SuprimeCam 
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24.7 


UltraVISTA 
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24.5 


UltraVISTA 


Y + J 


24.9 


UltraVISTA 


H 


24.0 


UltraVISTA 


K s 


23.8 


UltraVISTA 


3.6/im 


24.2 


Spitzer/IRAC 


4.5/^m 


23.8 


Spitzer/IRAC 



J imaging from WFCAM on UKIRT to report three potential 
z > 7 galaxies in COSMOS at surprisingly bright magnitudes 
(J < 23.3), some with proposed spectroscopic confirmation. How- 
ever, the power of the new UltraVISTA imaging is well demon- 
strated by the ease with which we can now show that all of these 
galaxies in fact lie at much lower redshift (see below). The main 
aim of this paper is, however, not only to check previous claims 
of z > 6.5 LBGs in the COSMOS field but more importantly to 
demonstrate that UltraVISTA, even in its first few months of data 
taking, has delivered images of the depth and quality necessary to 
produce a small but robust sample of luminous z ~ 7 galaxy can- 
didates, and to present this sample to the community for spectro- 
scopic follow-up. 

It is important to revisit why the discovery and study of 
rare luminous galaxies at these early epochs is of interest. First, 
while the study of faint galaxies may appear of more impor- 
tance for understanding reionisation, the errors on the faint-end 
slope (a) of the Schechter function parameterisation (4>(L) = 



d>*{L/L*y 



-L/L* 



) can be substantially reduced by a better de- 



termination of the form of the brighter end of the galaxy LF, due to 
the resulting improved constraints on the other key parameters (i.e. 
density normalisation <f>*, and characteristic luminosity L*). Sec- 
ond, it remains possible that, at some epoch, the exponential cutoff 
at the bright end of of the LF (still apparently present at z ~ 5 and, 
arguably, at z ~ 6; |McLure et al.|2009| |Willott et al.|2012| > may 
disappear and, for example, be replaced by a shallower power-law 
due to early inefficiency (or absence) of the physical processes (e.g. 
AGN feedback) purported to be responsible for the relative ineffi- 



ciency of star-formation in higher-mass dark-matter halos (e.g. Fin- 
|lator et al.|20"TTl but see also Trenti et al. 2010 and Jaacks et al. 
|2012[ >. Finally, quite apart from the issue of the form of the LF, the 
discovery of luminous z ~ 7 galaxies is of intrinsic interest, as it 
is often the most luminous galaxies which place the most stringent 
demands on models of galaxy formation (e.g. Benson et al. 2003) 
and z ~ 7 galaxies as bright as Y ~ 25 mag are clearly very at- 
tractive targets for follow-up near-infrared spectroscopy. 

This paper is structured as follows. In Section [2] we present 
the new UltraVISTA survey near-infrared data, and summarise the 



crucial supporting datasets: these comprise deep optical imaging 
(including the final CFHT Legacy Survey data, the long-public 
HST ACS i-band imaging, new ultra-deep Subaru z'-band imag- 
ing) and the mid-infrared (3.6 fim and 4.5 fim) Spitzer IRAC imag- 
ing obtained via the S-COSMOS survey l |Sanders et al.|2007) . In 
Section[3]we then describe the creation and subsequent progressive 
refinement of our galaxy sample based primarily on spectral en- 
ergy distribution (SED) fitting, and also show colour-colour plots 
to clarify the validity of our selection method and help expose the 
possible contaminant populations. Our final ten candidate z > 6.5 
galaxies are presented in Section [4] where we also discuss poten- 
tial contamination, with special emphasis on an extremely careful 
elimination of cool brown-dwarf (L, M, T) stars; as discussed by 
many authors (e.g. [SFanway et al. 2008 1, contamination by cool 
dwarf stars is a much more serious issue for ground-based searches 
for z ~ 6 — 7 galaxies than for ultra-deep HST WFC3/IR imaging 
surveys, both because extreme-redshift galaxies are generally un- 
resolved in even good-seeing ground-based imaging, and because 
the relative number of brown-dwarf stars to genuine high-redshift 
galaxies is much larger at brighter magnitudes (e.g. with a typical 
absolute magnitude J ~ 19, a T-dwarf star with an apparent mag- 
nitude of J = 24 still lies well within the galactic disc; see |Dunlop| 
|2012| >. Next, to place our results in context, in Section [5] we pro- 
vide a discussion/re-analysis of the previous claims of bright z > 7 
galaxies recently advanced by Capak et al. (201 1 ) and Hsieh et al. 
( |2012} , as well as the claimed highest redshift X-ray source can- 
didate presented by Salvato et aL~| j201l} . Then, in Section [6] we 
analyse the properties displayed by the stacked photometry of our 
top four galaxy candidates, and provide a very brief discussion of 
the implications of our results for the form of the LF at z — 7 (a 
detailed reanalysis of the z ~ 7 luminosity function is deferred 
to a future paper, where we aim to also incorporate data from the 
UKIDSS UDS survey). We conclude with a summary of our main 
results in Section 7. 

All magnitudes quoted are in the AB system (Oke & Gunn| 
|1983[ >, including data points from the dwarf-star literature, where 
Vega magnitudes and AB magnitudes are, unfortunately, frequently 
mixed. For the calculation of physical quantities we assume a 
ACDM cosmology with H = 70kms" 1 Mpc~ 1 , n m = 0.30 and 
Q A = 0.70. 



2 DATA 

The analysis presented in this paper is based on the first year of 
near-infrared imaging obtained by the UltraVISTA survey, in com- 
bination with the deep multi-band optical imaging obtained as part 
of the CFHT Legacy Survey (CFHTLS), new deep z'-band imag- 
ing obtained with the refurbished Suprime-Cam on the Subaru tele- 
scope, and other publicly-available HST and Spitzer data obtained 
as part of the Cosmological Evolution Survey (COSMOS) ( |Scov-| 
ille et al. 2007b). The production of multi-wavelength catalogues 
required the use of images on a common pixel grid and so the 
region of overlapping data defined by the CFHTLS imaging (de- 
scribed below in Section |2~2) was used to define the search area as 
shown in Fig.[T] In Table [T] we give the 5<r depths which we have 
calculated for each of the relevant imaging datasets (see caption for 
details). All images were matched to the astrometric grid of the Ul- 
traVISTA Y-band image using the IRAF package CCMAP. Finally, 
the images were resampled to match the pixel scale and image size 
of the CFHTLS D2 field using the IRAF package SREGISTER. 
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Figure 2. Colour-colour plots of our final ten candidates selected by SED fitting. The most secure high-redshift candidates are shown in blue. Candidates that 
are likely to be at high redshift, but where we cannot completely exclude the possibility they are dwarf stars or low-redshift galaxy contaminants, are shown in 
green. Finally, our our three least secure candidates are shown in red. The large uncertainty in the J — K s colour is a consequence of the relative shallowness 
of the K s imaging as compared to the Y+ J detection band (which also results in a few of our candidates being formally undetected in the K s image). The 
colours of M, L and T-dwarf stars are shown with a star symbol, and were obtained from photometric surveys I Knapp et al. 2004 , Burningham et al. 2010 1 and 
by calculating synthetic magnitudes from spectra obtained from several libraries as described in Findlay et al.H2012|. The grey tracks indicate the evolution 
of colours with redshift for high-redshift Lyman-break galaxies as synthesised from the theoretical SEDs of |Bruzual & Chariot (2003 ), assuming a constant 
star-formation history, and a range of other parameter values (Z = 0.2 Zq, Ay = 0.0 — 0.5, age = 50 — 500 Myr). Appropriate redshift steps are illustrated 



with grey circles situated along the tracks at z 
additional suppression of the flux blueward of A rcs t 



,6.5 and 7. Absorption by the Intergalactic Medium is applied with the Madau 1 1995) prescription, with 
1216 A at the very highest redshifts (parameterised as e~ 2 ^ 1+z ^^ from|Patel 2010 . 



2.1 Ultra VISTA near-infrared imaging 



2.2 CFHT optical imaging 



Over the next ~ 5 years it is planned that the Ultra VISTA sur- 
vey will provide deep near-infrared imaging covering the central 
region of the COSMOS field in four broad-band near-infrared fil- 
ters (Y, J, H, K s ) and one narrow-band filter (NB 1 18) to unprece- 
dented depths l McCracken et al. 2012'j^ . The survey was designed 
to commence with a "deep" programme (212 hours) providing 
fully-sampled imaging over a contiguous 1.5 deg 2 field, with the 
majority (1408 hours) of the observing time being subsequently 
devoted to the "ultradeep" programme, comprising substantially 
deeper imaging over four strips (0.73 deg 2 ) within the field. The 
first phase of "deep" imaging was completed by July 201 1, and re- 
leased to the public in fully-reduced form via ESO in Feb 2012. It 
is this year-one imaging that is utilised here, and the relevant cov- 
erage map and photometric depths are summarised in Fig. 1. The 
Ultra VISTA images use the COSMOS CFHT i*-band image from 
2003-2004 as the astrometric reference (McCracken et al.||20"T2"| 
Capak et al. 2007). For the present purpose of searching for high- 
redshift Lyman-break galaxies, an inverse-variance weighted stack 
of the Y-band and J-band Ultra VISTA images was created with 
the aim of increasing the sensitivity of the survey for objects with 
a near-flat (in /„) near-infrared spectral slope. 



Deep optical data are essential for the identification of genuine 
high-redshift galaxies, as it is vital to confirm that no flux is de- 
tected at wavelengths shortward of the putative Lyman-break. The 
COSMOS field benefits from extensive multi-band optical imag- 
ing obtained with both the CFHT and the Subaru telescope. In this 
study we use the deepest available data from the CFHT Legacy 
Survey, which provides imaging over a 1 deg 2 subsection of the 
COSMOS field centred on RA 10 h 00 m 28.00 s , Dec +2°12'30"; 
this defines the area of the survey appropriate for a consistent 
high-quality search for high-redshift optical drop-out galaxies. The 
CFHTLS T0006 releastj^Jin the deep field D2, provides data in the 
u*, g, r, ii, is ( "y" ) and z optical filters from the MegaCam instru- 
ment, with data in the it* -band obtained as part of the COSMOS 
survey (Capak et al.|2007] >. The ii and ii (or y-band) image dis- 
tinction is a consequence of a new i-band filter installed since Oc- 
tober 2007. However, since the filter transmission curves are similar 
and probe a region of the spectrum of z > 6.5 galaxies where we 
impose a non-detection condition, and the images are of compara- 
ble depth (to within 0.1 mag), we combined these images to form 
a deeper inverse-variance weighted stack, hereafter simply termed 
the i-band image. The astrometry of the CFHTLS T0006 data re- 
lease was based on the 2MASS catalogue rather than the standard 
COSMOS CFHT i*-band catalogue. Hence, the CFHTLS images 
were mapped onto the astrometric solution of the Y-band UltraV- 
ISTA imaging using the technique described in Section [2] with a 
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correction of typically 0.2", although the pixel scale was retained 
and used as a base for all other multi-wavelength images. 



2.3 HST/ACS I 8 i4-band imaging 

The COSMOS field has been imaged to single-orbit depth 
by HST/ACS in the Tgw-band through a HST Treasury Pro- 
gramme l|Scoville et al.||2~007a| |Koekemoer et al. 2007) |Massey| 
et al.|2010V Due to the large size of the 1.8 deg high-resolution 
mosaic (0.03"/pix). individual postage stamps of galaxy candi- 



and visually in- 
The /si-i-band 



dates were retrieved from the NASA/IPAC archivt 
spected (after smoothing) as described in Section 
image uses the same base astrometric reference as the COSMOS 
CFHT i*-band image used as the standard reference for subsequent 
imaging within the COSMOS field and the UltraVISTA data as de- 
scribed above, hence no further transformations were applied. 



2.4 Subaru Suprime-Cam z'-band imaging 

Deep z'-band imaging is crucial for the selection of galaxies at 
z > 6.5. Since 2009 we have exploited the new red sensitivity 
provided by the refurbished Suprime-Cam instrument on Subaru 
to obtain very deep z'-band imaging over the central 1 deg 2 of 
the COSMOS field (i.e. matching, as near as possible, the CFHT 
imaging). The imaging consists of four Suprime-Cam pointings, 
each with > 15 hours of exposure time, rotated by 90 deg with 
respect to each other. This yields a final .z'-band mosaic reaching 
a minimum 5a depth of z = 26.3 (in a 2-arcsec diameter aper- 
ture), with the deepest panel reaching 0.3 mag deeper. To create the 
mosaic, the astrometry of each of the four pointings was matched 
to that of the Y-band UltraVISTA imaging using the IRAF pack- 
age CCMAP, background-subtracted using SEXTRACTOR ( |Bertin &| 
|Arnouts|1996| >, and the zero points were equalised to take into ac- 
count the different exposure time of each pointing. A mosaic in the 
native pixel size (0.202"/pix) was produced using SWARP (Bertin 



|et al.|20 02), where overlapping sections were combined with the 
WEIGHTED keyword using weights taken from an rms map created 
by SEXTRACTOR. A science image for use with the UltraVISTA 
datasets, overlapping the CFHT Legacy Survey image region with 
a pixel scale of 0. 186"/pix, was created using the IRAF package 
SREGISTER. The full Subaru mosaic overlaps the CFHTLS field 
that forms the search area in this paper, with some sections extend- 
ing beyond as can be seen in Fig.[T] 



2.5 Spitzer IRAC mid-infrared imaging 

The full 2 deg 2 COSMOS field is covered by publicly-available 
mid-infrared data from Spitzer IRAC and MIPS obtained as part 
of the S-COSMOS survey {Sanders et al.|2007| >. All channels were 
imaged with an integration time of 1200 seconds. The IRAC 5.8/im 
and 8.0/im data are too shallow (rriAB = 21.3 and 21.0 respec- 
tively, Sanders et al. 20071 for effective use in the present study, 
beyond visual inspection of our final candidate sample (see Sec- 
tion [3j. However, the IRAC 3.6 fim and 4.5 /im images reach 5a 
depths in a 2.8"diameter aperture of 24.2 and 23.8 mag respec- 
tively, are sufficiently well-matched to the UltraVISTA H and K s 



4 |http : //irsa ■ ipac ■ caltech . edu/data/COSMOS/index_ | 
|cutouts . html| 



depths to be of potential use in the galaxy selection process. In par- 
ticular, low-redshift galaxy contaminants are predicted to be par- 
ticularly bright at these wavelengths, in contrast to the much flatter 
near-infrared slope anticipated for genuine high-redshift galaxies. 
The Spitzer IRAC photometry was obtained in a 2.8"diameter aper- 
ture, from images that had been background subtracted with a large 
filter length using the GLOBAL keyword in SEXTRACTOR and then 
matched to the common pixel grid and area defined above. In ac- 
knowledgement of the limitations of photometry from the confused 
IRAC imaging, the 3.6 and 4.5 /im photometric errors used in 
the photometric redshift fitting described below were set to a mini- 
mum level of 20% of the measured flux density. 



3 CANDIDATE SELECTION 

3.1 Initial detection, photometry and depth analysis 

The primary catalogue was created using SEXTRACTOR 
v2.8.6 {Berlin & Arnouts||1996| > on the Y+J stacked image, 
with photometry in the additional filters collected in the dual- 
image mode. Objects were also added to the catalogue from Y- and 
J-band selected catalogues, to ensure no galaxies were missed. 
However, all objects retained in the final sample were in fact 
selected from the Y+ J catalogue. All optical and near-infrared 
photometry was measured in a 2"diameter circular aperture, which 
corresponds to 70% enclosed flux (for a point source) in the 
Y-band imaging. 

The global limiting magnitude for each image was calculated 
by finding the sigma-clipped standard deviation of the flux enclosed 
within 2"(or 1 .2" where relevant) diameter circular apertures placed 
on the image in a grid. Apertures that contain flux from an ob- 
ject were removed using the segmentation-map output produced by 
SEXTRACTOR, where aggressive parameter values were chosen to 
ensure detection of the majority of low-significance objects (ex- 
ceeding 2a significance in the detection procedure). For the final 
sample of objects presented in this paper, we undertook a local es- 
timate of the photometric errors by finding the standard deviation of 
the 50 closest blank apertures placed around the object, with 2.5er 
clipping. 

The 5a depth of the Y+ J image is 24.9 mag, and so, for sim- 
plicity, we cut our catalogue at Y+ J < 25.0, resulting in an initial 
near-infrared selected sample of 175075 sources. To confine our 
search to potential objects at redshifts z > 6.0, we then applied 
a non-detection (< 2 a) criterion in all optical filters blueward of 
and including the i-band. Removing artefacts around bright stars 
and the small region of missing UltraVISTA data in the corner of 
our central survey area, resulted in a catalogue of 146 objects with 
Y+ J < 25.0 and non-detections in all optical bands except for the 
CFHT z-band and/or Subaru z'-band imaging. 



3.2 Sample refinement via SED fitting 

In contrast to other studies which have endeavoured to select bright 
z ~ 7 galaxies from ground-based data, we do not attempt to define 
areas in colour-colour space to isolate high-redshift galaxies. Fig. [2] 
clearly illustrates the difficulties in separating high-redshift galax- 
ies from dwarf-star contaminants based on z', J and K s photome- 
try alone. The addition of Y-band data is crucial (as shown in Fig. 
2b), but again the use of only two colours does not make optimal 
use of the available multi-wavelength data. Therefore, as in our pre- 
vious studies (e.g. McLure et al. 2009; 2010; 201 1) we have used a 
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Figure 3. The colours of our high-redshift galaxy candidates (shown as 
the coloured circles) using the Spitzer IRAC information, where the most 
secure candidates are shown in blue, the robust/contaminant category are 
shown in green and the insecure category are shown in red. The colours of 
M, L and T-dwarf stars are shown with grey star symbols from Kirkpatrick 
|et al.|(2~011) and |Patten et al.| l2006l, where the error bars are all below 
0.05 mag. Predicted colours of low-redshift contaminant galaxies (1.3 < 
z < 1.9) are shown by the hatched area, where reddening increases to 
redder J— 3.6/xm colours (from A\r — to ~ 2 at the right-hand edge of the 
plot). The thin navy wedge represents the predicted colours of high-redshift 
galaxies using the models and parameter ranges as in Fig. [5] Flux densities 
below the 2cr level were set to the 2a depth and the resulting colour limits 
plotted as an arrow. Note that T-dwarf stars have the reddest 3.6itm - 4.5/im 
colours, whereas L and M types all have slightly bluer colours and occupy 
the lower horizontal band of stars. 



full SED-fitting analysis to derive the redshift-probability distribu- 
tion for each galaxy candidate. We employed the Bruz uaT& Chariot] 
(2003) stellar evolution models, considering models with metallic - 
ities ranging from solar (Zq) to l/50th solar (0.02Zq). The star- 
formation histories considered were instantaneous bursts, constant 
and exponentially declining with characteristic timescales in the 
range 50 Myr < r < 10 Gyr. The ages of the stellar populations 
models were allowed to range from 10 Myr to 13.7 Gyr, but were 
required to be less than the age of the Universe at each redshift. 
Dust reddening was described by the Calzetti et al. (20001 attenu- 
ation law, and allowed to vary within the range 0.0 ^ Av ^ 4.0. 
Inter-galactic medium absorption short-ward of Lya was described 
by the |Madau| ( |1995fr prescription, and a |Chabrier| ( |2003^ IMF was 
assumed in all cases. Further details can be found in lMcLure et al.l 
< |201 1| >. Specifically for this study we have included detailed fit- 
ting of the spectra of L, M and T dwarf stars, where the reference 
spectra from the speX librar^ were used for each spectral type 
from M4 through to T8. Only for the final sample of galaxies pre- 
sented in Section [4] did we introduce the equivalent width (EW) 
of the Lya -line in the SED-fitting analysis as an extra free pa- 
rameter, with the rest-frame equivalent width within the range A 
< EW < 240 A. 

SED fitting was performed using all the multi-wavelength 
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data, including the CFHTLS z-band data (although it is signifi- 
cantly shallower than the Subaru z'-band image, it has a slightly 
bluer effective wavelength and thus adds some extra spectral res- 
olution to the fitting process). The IRAC photometry was only in- 
cluded in the final stages of candidate selection due to the large 
uncertainties involved and the high number of objects where the 
photometry is confused. All near- and far-infrared photometry was 
corrected to the 84% enclosed flux level of the z'-band imaging 
for the SED fitting, where the enclosed flux within a 2"diameter 
aperture was calculated from the point spread function obtained by 
taking the median of unsaturated stars extracted from each image. 

The output redshift-x 2 distribution from the SED fitting (with- 
out the IRAC photometry) was used to determine which of the 
146 candidates were consistent with being at high redshift. We ap- 
plied the criteria that the object must have an acceptable solution 
(\ 2 < 11.3, which corresponds to 2a significance given 5 degrees 
of freedom) at z > 6.0, and that the z > 6.0 template must be pre- 
ferred (Ax 2 > 0.0) over the alternative low-redshift galaxy solu- 
tion (usually produced at a redshift where the putative Lyman-break 
can be interpreted as a Balmer or 4000A break). These conditions 
resulted in a reduced sample of 68 objects consistent with being at 
z > 6.0. Further inspection of the i-band image, in combination 
with the HST/ACS 7si4 image and an inverse-variance weighted 
stack of the optical bands up to and including i, resulted in the re- 
moval of a further 21 objects with weak detections in these bands 
implying z < 6.5. 

As a result of the SED fitting and manual optical checks, the 
sample of potential z > 6.5 galaxies was reduced to 47 objects. 
However many of these still had formally-acceptable solutions at 
much lower redshift. At this point the IRAC 3.6 /im and 4.5 pirn 
measurements were incorporated into the SED-fitting analysis, in 
effect to remove the red dusty low-redshift galaxy interlopers which 
display a much redder J — 3.6/^m colour than genuine high-redshift 
galaxies. In practice the IRAC photometry was simply incorporated 
into the full SED fitting, but the process is illustrated in colour- 
colour space in Fig. [3] which shows the (hashed) area at J — 3.6 pirn 
> 1 where low-redshift dusty interlopers are generally found. At 
this stage all the observed optical-to-near-infrared SEDs were also 
fitted to the SEDs of L, M, T dwarf stars and objects were rejected 
if the dwarf-star solution was formally preferred over the high- 
redshift galaxy solution. The lack of mid-infrared spectroscopy of 
dwarf stars prevented the simple incorporation of the IRAC mea- 
surements into the dwarf-star SED-fitting process, but we were able 
to confirm a number of ambiguous dwarf-star solutions on the ba- 
sis of the position of the object on the colour-colour plots shown 
in Fig.|2]and Fig. [3] The final result of this IRAC-based cull was a 
remaining sample of 17 potential z > 6.5 galaxies. 

The final stage of sample refinement involved detailed re- 
investigation of every single galaxy image, and repetition of the 
dwarf-star fitting using photometry taken in a smaller 1.2"diameter 
circular aperture (with the expectation that, in the case of a real 
star, the fit should improve with the smaller aperture photometry). 
Galaxy colours were also calculated from smaller-aperture pho- 
tometry for comparison. As a consequence of these final checks, 
three objects were excluded because of low-level flux detected in 
the optical stack, in combination with insecure photometry that re- 
sulted in the objects drifting substantially around the colour-colour 
diagrams when smaller apertures were used. A further two objects 
were only detected in the V-band image, and closer inspection re- 
vealed them to be part of a faint diffraction halo around a bright 
star. Finally two objects were confirmed as secure z > 6 galaxies, 
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but their redshift-probability distributions indicated that they could 
not lie above z = 6.5. 

The final outcome is a sample of 10 objects which have sur- 
vived as credible candidate Lyman-break galaxies at z > 6.5, and 
which we deem worthy of presentation in this paper. The magni- 
tude errors and SED fits for these objects presented in subsequent 
sections are all based on the local error estimates (although these 
are all formally consistent with the global error estimates). We vi- 
sually inspected the 5.8 /im and 8.0 /im SpitzerfiRAC imaging for 
our final 10 candidates, finding no evidence for even low-level flux 
(as expected given the magnitudes predicted by the SED fits, and 
the depths of the S-COSMOS IRAC imaging at these longer wave- 
lengths; see Section |23) . 



4 CANDIDATE z > 6.5 GALAXIES 

Table|2]details the photometry of the final ten z > 6.5 galaxy can- 
didates, where we have grouped them into three categories depend- 
ing on the statistical security of their z > 6.5 solutions. The final 
galaxy and dwarf-star SED fits for these ten candidates are pre- 
sented in Fig. [5] Their photometric redshifts and derived physical 
properties such as stellar masses and star-formation rates are pre- 
sented in Table [3] 

4.1 Category 1 - Robust 

Our four most secure z > 6.5 galaxy candidates have completely 
unacceptable alternative fits for either a low-redshift dusty galaxy 
or a late-type dwarf star from the SED fitting analysis, and are all 
very well described by a high-redshift galaxy template at z > 6.5, 
as can been seen in Fig. [5] 

Object 277912 has z — Y = 2.5 and a flat near-infrared SED 
through to the IRAC bands, with a best-fitting spectral template 
corresponding to a galaxy with z p hot = 6.971q'q7. As can be seen 
from the postage stamps presented in Fig. [4] a low-redshift galaxy 
is present ~2.5" away from the position of our high-redshift can- 
didate in the Y-band imaging. The companion has a photometric 
redshift of z phot = 0.94 ± 0.10 from the COSMOS Photometric 
Redshift Catalogue filbert et al.|2008) and z phot = 1.18 ± 0.10 
from our own fitting of the 2 diameter aperture photometry. The 
X 2 distribution for object 277912 indicates that the low-redshift so- 
lution is not only extremely unlikely, but would also lie at z ~ 1.7. 
Hence it seems clear that this object is at high redshift and not as- 
sociated with the low-redshift foreground galaxy. Unfortunately the 
presence of the foreground galaxy means that the IRAC photome- 
try, where the FWHM of the PSF is ~ 1.7", is confused, and there- 
fore the stellar mass estimate is uncertain. Allowance for a possible 
contribution from Lyman-a emission in the SED fitting does not 
change the photometric redshift, with zero Lyman-a emission re- 
maining the preferred option. 

Object 155880 also has a blue spectral slope in the near- 
infrared bands and is clearly detected in the deep Subaru z'-band 
imaging giving a z p hot = 6.781q'o9- We stress that the z'-band de- 
tection (z' — 26.1 ± 0.1) is achieved because the Subaru Suprime- 
CAM z'-band imaging is so deep and the filter extends to ~ lfim; 
the measured spectral break is still large z' — Y > 1.5, although 
not large enough to warrant inclusion as a high-redshift candi- 
date in some traditional colour-colour cut selection methods (e.g. 
z' — Y/J > 2). Fitting with the addition of a Lya line results in 
a slightly higher-redshift fit with z p hot = 6.98 (with a Lya-line 
EWq = 40 A). The best-fitting stellar template of an earlier type 



T-dwarf star (T3) can produce the z — Y colour, but not the other 
near-infrared colours, with the result that the high-redshift galaxy 
solution is the only acceptable SED fit. The IRAC 3.6/im - 4.5/xm = 
0.6 colour is redder than would be predicted directly from constant 
SFR models, but there is evidence (Curtis-Lake et al, in prepara- 
tion) that a multiple-component star-formation history or a signifi- 
cant contribution from nebular emission can reproduce this spectral 
shape. Inspection of the smoothed ACS TgM-band imaging shown 
in Fig. [4] shows a low-significance (2.3<r) detection 0.5" from our 
object. However, given that our astrometric accuracy is better than 
0.2", we conclude that this detection is unrelated to our candidate. 

Object 218467 has the bluest Y — J colour of all our candi- 
dates as shown in the colour-colour plot in Fig. [2] which results 
in unacceptable SED fits for both the low-redshift and stellar tem- 
plates, as neither can reproduce both the magnitude drop between 
z and Y and the blue spectral slope at longer wavelengths. The J- 
band detection is below the 5 a limit of the data and hence there is a 
large uncertainty in the extremely blue Y — J = —0.5 colour (bluer 
than predicted by our range of models; see Fig.|2j. The candidate is 
formally undetected in the current K s and //-band imaging which 
leads to a large uncertainty on these data points in the SED fitting. 
However, when a smaller 1.2" diameter aperture is used, marginal 
detections produce a flat spectral shape as seen in the dwarf-star 
SED fit figure. A faint detection in the z'-band is consistent with the 
z p hot = 7.04^Q'Jg within the errors, although the detection here is 
best fitted with inclusion of a strong Lya line of EWo = 110 A at 
the very red end of the Subaru z'-band filter at z p i lot = 7.20. Fig. [4] 
shows a detection in the IRAC 3.6 fim image consistent with a flat 
SED extending from the near-infrared data, but the 4.5 image 
appears to have a negative hole at the position of our object. The 
implied stellar mass of 218467 is the lowest of the four robust can- 
didates, indicating that it may be underestimated as a result of the 
potentially spurious lack of 4.5 fim flux. Alternatively, the implied 
blue 3.6^im - 4.5/im slope could indicate the presence of strong 
nebular emission lines. 

Finally, object 61432 has a flat near-infrared spectral slope 
and has a strong spectral break, with only a marginal detection in 
the z'-band. The somewhat red 3.6pim - 4.5/im colour might ap- 
pear difficult to reconcile with a high-redshift SED as illustrated in 
Fig. [3] but the error bars are large and neither a late-type dwarf star 
nor a low-redshift galaxy provide remotely acceptable solutions. 
As with the first object above (object 277912) allowance for a pos- 
sible contribution from Lyman-a emission in the SED fitting does 
not change the photometric redshift, with zero Lyman-a emission 
remaining the preferred option. 

4.2 Category 2 - Robust/Contaminant 

The three candidates that make up category 2 are all still consis- 
tent with being at high redshift, and the z > 6.5 galaxy solution 
is still formally preferred. However, with the current data we can- 
not exclude the possibility that these objects could be either at low 
redshift or galactic dwarf stars. 

The best-fitting SED for object 277880 is a galaxy at z p hot = 
6.77+Q gg. However, because of a strong z'-band detection in 
combination with a relatively red spectral-slope through the near- 
infrared bands, the low-redshift solution at z = 2.0 cannot be com- 
pletely excluded. A stellar fit of type T3, although not as good 
as that of the high-redshift galaxy SED fit, cannot be completely 
excluded either. However, the flat Y — J colour is hard to repro- 
duce with either a low-redshift galaxy or dwarf-star fit, and so the 
z > 6.5 galaxy solution is still favoured for 277880. Allowance 
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Figure 4. Multi-band postage-stamp images of the four most secure z ~ 7 galaxies presented in this paper, with the detection image (Y + J) shown on the 
left. A stack of these four objects is presented in the bottom row. The stamps are 10" X 10", with North up, and East to the left. Each filter stamp has a linear 
grayscale, with the lower level (white) set to 2<r below the background level and objects (black) saturated at 4.5<r above the background level. To aid the visual 
detection of low-level flux indicative of a z < 6.5 contaminant object, the 7gi4 stamp has been smoothed using the IRAF GAUSS package, with a standard 
deviation of 5 pixels. 



Table 2. Photometry of the ten z > 6.5 galaxy candidates, grouped into the three categories described in the main text. From top to bottom the top four 
galaxies are classed as "Robust" (followed by the photometry for a stack of this subset of galaxies), the next three as "Robust/Contaminant", and the final 
three can be regarded as "Insecure". The optical and near-infrared magnitudes were measured within a 2"diameter circular aperture and corrected to 84% 
enclosed flux, and the IRAC magnitudes were measured within a 2.8"diameter circular aperture and corrected to 84% enclosed flux. As a consequence of the 
correction up to 84% enclosed flux, the Y- and J-band magnitudes presented here appear typically ~ 0.2 brighter than the quoted 2-arcsec diameter 5cr limit 
(see Table[TJ. Wherever the measured flux density lies below the 2a limit in that band, the magnitude has been replaced here by the 2cr limiting magnitude 
depth of the appropriate image, and is given as a lower limit on the apparent magnitude. 
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Figure 5. SED fits for each member of the final sample of ten high-redshift galaxies presented in this paper. Where the candidate has a non-detection (below 
the lcr level) in a given band, that point has been plotted with an arrow where the tip represents the lcr limiting depth. The best-fitting high-redshift galaxy 
template (without Lya) is shown in blue in the left hand panel, with the best-fitting alternative low-redshift template shown in red. The inset shows the 
redshift-x 2 distribution produced in the fitting process, where the redshift range 6.5 < z < 7.5 targeted in this study is highlighted in grey. The best-fitting 
low-redshift galaxy templates have a redshift around z ~ 1.2 — 1.6, due to the fact that a 4000A/Balmer break in the model galaxy spectra can sometimes 
reproduce the observed spectral break at ~ lfj,m if the data have insufficient dynamic range. For each object the best-fitting stellar template is shown in the 
right-hand panel, where photometry measured in a 1.2"diameter aperture has been used. The inset for the dwarf-star plots shows the \ 2 distribution versus 
stellar type, with effective temperature decreasing from M4 through to T8. The resulting best-fit photometric redshift, stellar type and derived parameters are 
detailed in Table[5] 
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Figure 5. Continued. 
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Table 3. Results from SED fitting of the ten z > 6.5 galaxy candidates. The photometric redshift for the high-redshift solution is given, along with the raw 
Xg a p f° r fitted galaxy templates with or without a contribution from Lya line emission. The dust-extinction (Ay) and metallicity (Z) values given in columns 
4 and 5 are taken from the best-fitting model across all star-formation histories considered (but without Lya ). For the fits involving Lya we give simply the 
inferred photometric redshift and the rest-frame equivalent width (EWo) of the Lya line as selected by the best-fitting model. By definition, the X 2 gsi values 
for the Lya fits are either the same as, or lower than the x ga i values given for the continuum-only fits as listed in column 3, while the photometric redshifts 
inferred from the Lya fits are (inevitably) either the same as, or greater than those given in column 2. Note t hat a s a result of including the IRAC photometry 
in the final stages of SED fitting, the x 2 values here can exceed the x 2 < H-3 condition stated in Section [3~2| Columns 9 and 10 give the key information 
for the best-fitting dwarf-star SED fits (note that, because of the different number of degrees of freedom, the values of Xgtcllar cannot b e directly compared 
with the corresponding values for the galaxy fits). Finally, in the last 5 columns we give derived physical parameters for these objects without Lya in the 
fitting. For stellar mass (M*) we simply quote the range of possible stellar masses as calculated from tau, constant or burst star-formation models, while 
the tabulated star formation rate (SFR) and specific SFR (sSFR) are calculated from the best-fitting constant star-formation model. The Ay value for the 
best-fitting constant star-formation model is also presented, where the most extreme SFRs are inevitably coupled with the larger dust-reddening values (see the 
discussio n in Section[6). The final colum n gives the SFR estimated directly from the (observed) rest-frame UV luminosity (i.e. assuming zero dust extinction), 
using the Madau, Pozzetti, & Dickinson 1 1998 1 prescription (dividing by a factor of 1.8 to convert from a Salpeter to a Chabrier IMF). 
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for a possible contribution from Lyman-a emission in the SED fit- 
ting does not change the photometric redshift, with zero Lyman-a 
emission remaining the preferred option. 

For object 268511, the SED fitting again indicates that the 
high-redshift solution is the best fit, with the strong z 1 — Y > 2.6 
break resulting in z p hot = 7.071q '\q and excluding all stellar fits 
except for the reddest T-dwarfs. Inspection of the postage-stamps 
for this object, however, shows that it is only clearly seen by eye in 
the y-band where it is compact, and therefore it could be a transient 
object or a T8 star. In this case, allowance for a possible contribu- 
tion from Lyman-a emission in the SED fitting raises the photomet- 
ric redshift to z p hot = 7.22, with an implied Lyman-a rest-frame 
equivalent width of EWo = 50 A. 

Finally, 271105 has a best fitting high-redshift SED with 
•Zphot = 6.621q'^, with a higher-redshift solution (z phot = 6.97) 
possible with the introduction of Lya emission with EWo = 70 A. 
A clear detection in the z'-band does, however, allow the stellar 
fit for a T3 dwarf to recreate the spectral break. However, the red 
J— 3.6/im colour of 271 105 is inconsistent with a T-dwarf star (the 
object is well detected and unconfused in all near-infrared bands). 
Hence, although we cannot rule a stellar solution, it still seems 
highly likely that 271 105 is indeed a high-redshift galaxy. 



4.3 Category 3 - Insecure 

The final three candidates all have good SED fits to a high-redshift 
galaxy template, with photometric redshifts z < 7.0 due to clear z - 
band detections. However, for these objects an M-dwarf stellar fit 
cannot be completely excluded, and some of the photometric red- 



shift estimates suggest a redshift just below z = 6.5. For both these 
reasons we classify these objects as insecure z > 6.5 galaxies. 

Object 95661 has a high-redshift galaxy fit with z p hot = 
6.13^0 27' an d so still has some probability of being at z > 6.5 
(although here the introduction of Lya emission does not increase 
the inferred redshift). A dwarf star of type M6 can also reproduce 
the flat near-infrared photometry as well as the spectral break, and 
so cannot be excluded. 

Object 28400 has a high-redshift galaxy fit with z p hot = 
6.34^019 but in this case is actually better fitted with the inclusion 
of a Lya emission line with EWo = 70 A increasing the inferred 
redshift above the z — 6.5 threshold to z p hot = 6.65. The IRAC 
photometry shows a bright detection in the 3.6 /im image but noth- 
ing in the 4.5 /im data despite the similar depth, implying a break 
of at least 3.6/im - 4.5/xm < —1.4 (taking the 2a limit for 4.5 (im). 
A possible explanation would be an M-dwarf star which can pro- 
duce this blue colour within the error bars as shown in Fig. [3] One 
alternative explanation for the blue 3.6/im - 4.5/im colour would 
be nebular emission. 

Finally, object 2233 has a clear z'-band detection, but is weak 
in Y+ J and the other infrared bands. The high-redshift galaxy fit 
is comparable in quality to that of a M6 star and the weak near- 
infrared photometry allows the inclusion of a Lya line of EWo = 
150 A which raises the photometric redshift to Zphot = 7.03. 
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5 COMPARISON WITH PREVIOUS STUDIES 

Several other ground-based studies have attempted to determine, or 
at least place a meaningful limit on the very bright end of the galaxy 
UV luminosity function at z ~ 7, using deep near-infrared photom- 
etry and colour cuts to attempt to exclude low-redshift galaxies and 
cool galactic dwarf-star contaminants. 



5.1 |Capak et al.|p011) 

Of particular relevance to this paper is previous work within 
the COSMOS field undertaken by |Capak et ai] poll) . They 
utilised the deepest multi-wavelength data over the full COS- 
MOS 2 deg 2 field available at the time, in particular optical data 
from Subaru including the i + and z + -bands, near-infrared data 
in the J-band from UKIRTAVFCAM and H, K s imaging from 
CFHT/WIRCam < |Bielby et alT][20TT] |McCracken et~aLl[20T0| ). 
Galaxy candidates were selected in the K a -ban& with a 5a cut 
in the J and K s bands of 23.7 mag (3"diameter aperture) applied 
along with the colour conditions z' — J ^ 1.5, K s — 4.5/im > 
and J — K s > 0. They found three objects consistent with be- 
ing at z > 7 as defined by their colour-colour criterion, and pre- 
sented follow-up spectroscopy which, along with SED fitting, re- 
sulted in the rejection of one object (C2 below) as most likely lying 
at much lower redshift. Since the near-infrared data from Ultra V- 
ISTA extends to around one magnitude deeper than that exploited 
by |Capak et al.| {201 1[ > in the J-band (with improved seeing) and 
also provides F-band imaging that greatly improves the separa- 
tion of contaminant populations (as described previously), we have 
extracted new photometry for all three of these objects and under- 
taken SED fitting in the same manner as described above for our 
own, new deeper sample of objects. Since none of the Cap ak et"aI7| 
( |20 1 1| > candidates actually made it into the sample described above, 
it should not be a surprise that, as described below, we fail to con- 
firm that they lie at z > 6.5. 

The first object presented by Capak et al. ( 201 1 ), hereafter CI, 
in fact lies just outside the area of overlapping deep data utilised 
in this paper and, as a consequence, the CFHTLS and new deeper 
Subaru z'-band imaging is unavailable. However, CI is still within 
the 1.5 deg 2 of the Ultra VISTA imagingthat extends beyond the 
central ~ 1 deg 2 area, as seen in Fig. Ill With the deeper near- 
infrared data, most importantly the F-band, we obtained photom- 
etry for this object to check the consistency of the photometric 
redshift calculated in Capak et al. (2011). Near-infrared photom- 
etry was measured within a 2"diameter aperture using SEXTRAC- 
TOR in dual-image mode with selection in the A" s -band, where 
the candidate was brightest. For the optical bands, cutouts from 
the COSMOS website in the Bj, g + , Vj, r + , i + and z + -bands 
were extracted and the magnitudes (measured using GAIA) and 
errors (taken from the 5a limits presented in Capak ^t al.||201l) 
were included in the SED fitting analysis. We find CI to be unde- 
tected in Y, and to have J = 24.4 ± 0.2, H = 23.3 ± 0.1 and 
K s — 23.2 ±0.1 (within a 2"aperture with a correction to 84% 
enclosed flux). These values can be compared to the Cap ak et al.| 
pOTT) magnitudes of J = 23.21 ± 0.05, H = 23.09 ± 0.12 and 
K 3 = 22.14 ± 0.02 within a 3"aperture. The large differences 
in photometry cannot be explained by aperture corrections, as the 
3 "aperture used by Capak et al. {201 1| was predicted to enclose 
80% of the flux in the PSF-matched J, H, K 3 imaging, implying 
the new magnitudes should be slightly brighter or comparable. An 
offset in photometry between the WIRCam and Ultra VISTA data 
has been identified in the UltraVISTA data release documentation 




Figure 6. SED fits for the three z > 7 galaxy candidates from |Capak| 
|et al.| {2011} using our revised photometry from the deeper UltraVISTA 
near-infrared imaging and the new 2' -band Subaru imaging. The UltraV- 
ISTA data extends to around one magnitude deeper in the J-band that the 
previous imaging available in the COSMOS field. In addition, it provides 
previously unavailable deep y-band imaging, invaluable for removing low- 
redshift galaxy contaminants and dwarf stars. The blue line shows the best 
fitting high-redshift model (with z > 4) and the red line shows the best 
low-redshift fit, where the photometric redshift and x 2 values for these 
two solutions are Capakl: z p hot — 

2.1(7.6), x 2 = 9.0(26.0), Capak2: 
Zpkot = 1-4(6.9), x 2 = 5.0(42.2) and Capak3: z phot = 3.4(12.6), 
X 2 = 2.2(9.1). We have extended the redshift range of the Capak3 fit 
above to display the high-redshift model solution at z ~ 13, which al- 
though being a formally acceptable fit has an implied star-formation rate 
and stellar mass that make it unrealistic at this extreme redshift. Hence, 
with our new photometry it can be seen that none of these objects has a 
robust high-redshift fit, but rather all lie in the redshift range z ~ 1.5 — 3.5. 
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(see also |McCracken et al.|2012) , implying that the UltraVISTA 
data is 0. 15 mag fainter than the previous imaging, which is thought 
to be a problem with the COSMOS data rather than with VISTA; 
however this effect still cannot explain the change in the observed 
near-infrared colours. The IRAC imaging of CI is highly-confused 
with a nearby object and so no attempt was made to measure this 
by either ourselves or byjCapak et al. (2011). The large compan- 
ion object is a low-redshift interloper at 2 p hot = 0.35 from the 
COSMOS catalogue, which may have affected the accuracy of the 
3"photometry. Our SED analysis of CI, with the fit shown in Fig. [6] 
yields no acceptable high-redshift solution, but rather reveals this 
object to be a red galaxy at z p hot — 2.1. 

The second object, C2, was already identified by CapakeFaT] 
< |20 1 1 ) > as a likely low-redshift interloper at z p hot = 1.59, which 
we here confirm with our SED fit to the UltraVISTA data shown 
in the second panel of Fig. [6] This galaxy candidate was selected 
in our own preliminary catalogue, but excluded at a relatively early 
stage as a low-redshift galaxy because of the clear i-band detection. 

Surprisingly, the third object, C3, was not selected in our anal- 
ysis despite the claim by |Capak et al.| l |201 1| > that it has J = 
23.1 ± 0.1. Visual inspection of the UltraVISTA J-band image of 
this object in fact reveals little, if any, evidence of detectable flux, 
however the object is real as it is clearly visible in the 7\~ s -band. Ac- 
cordingly, photometry from a K s -band selected catalogue for the 
object was collected and the SED fit performed which, as shown in 
the third panel of Fig.|6] reveals the galaxy to lie at z ~ 3.5. 

This analysis therefore highlights the power of the new deeper 
and better-quality UltraVISTA data over the previously available 
near-infrared data in the COSMOS field, and confirms that the new 
objects uncovered here from the Y+ J UltraVISTA imaging are the 
brightest credible z > 6.5 galaxy candidates in this area of sky. 

5.2 ISalvato et al.|p011) 

Another particularly relevant high-redshift candidate within the 
field is the object CID-2550 identified by |Salvato et ah] pOTT) , 
that they claim could potentially be the highest redshift X-ray- 
selected source found to date. By SED fitting to the optical and 
near-IR counterpart to a Chandra X-ray source (from the Chandra- 
COSMOS catalogue), they found a best fit photometric redshift of 
Zphot ~ 6.84. The previous magnitudes of this object CID-2550 
are i « 26.6, z = 25.4, J = 23.6, H = 23.8,^ = 23.0 with 
the object becoming brighter still in the Spitzer channels. However, 
despite the extremely bright J-band magnitude reported by |Salvato] 
|et al.|p01 1[ >, CID-2550 was not present in our final sample. CID- 
2550 is present in our original Y + J selected catalogue as you 
would expect, but was immediately excluded as it is unexpectedly 
below the 5<7-limit in the detection image (with Y + J = 26.9). In- 
spection of the UltraVISTA imaging and associated datasets shows 
the candidate to be clearly detected in the H and K s -bands, with 
faint detections in the i, z , Y and J-bands. As the object is bright- 
est in the K s -bmd we used A" s -band selected magnitudes and per- 
formed SED fitting as for our own candidates. Figure [7] shows 
our best fit to this object at z p hot = 2.6, where the decrease in 
the photometric redshift compared to the previous determination 
by |Salvato et al.| ( |2011[ l is a consequence of the reduced J-band 
magnitude and the inclusion of the y-band, along with reduced 
errors throughout the near-IR that show the spectrum to be gradu- 
ally rising rather than exhibiting a sudden break in the photometry. 
Given that the candidate is an X-ray source, we conclude that it is 
most likely a dusty active galactic nucleus at a lower redshift of 
%>hot = 2.6. 




Figure 7. SED fit for the candidate z ~ 6.84 X-ray-selected source 
from Salvato et al. 1 20 1 1 i using our revised photometry. The blue line shows 
the best fitting high-redshift model (z p hot = 6.7, x 2 = 22.9) and the red 
line shows the best low-redshift fit (z p hot = 2.6, x 2 = 1.9). Due to the 
improved near-IR photometry, which shows the candidate to have a gradu- 
ally rising SED rather than presenting a strong break, we conclude that the 
object is in-fact a lower redshift, dusty interloper. 



5.3 |Hsieh et al.U2012) 

|Hsieh et al.| ( |2012) recently reported a search for z$$o -dropout 
galaxies within the Extended Chandra Deep Field South (ECDFS), 
with objects detected in J-band data obtained as part of the Taiwan 
ECDFS Near-Infrared Survey (TENIS). Boxes in colour-colour 
space were used to select high-redshift galaxies, in particular the 
Spitzer IRAC channels were introduced in an attempt to elimi- 
nate low-redshift galaxy and dwarf-star contaminants. One can- 
didate high-redshift galaxy with J = 25.12 ± 0.23 (total mag- 
nitude), TENIS-ZD1, was selected from the 0.25 deg 2 surveyed, 
with a weighted photometric redshift of z p i lot = 7.822 jlo'°||. 
The very large errors on the derived photometric redshift are a 
consequence of the proposed Lyman-break occurring between two 
widely-separated filters, the zsso- and J-band, due to the unfor- 
tunate lack of Y-band imaging. It is this lack of y-band imaging 
which also makes it hard to exclude dwarf-star contaminants us- 
ing the techniques discussed above (Fig. |2^), leading Hsieh et al. 
( |2012} to rely on the IRAC colours in an attempt to achieve this. 

The infrared colours of dwarf star contaminants utilised 
by |Hsieh et al.H20T2| were calcu lated from the AMES-dusty dwarf 
star models ( |Allard et al.|2001| >. However, with the advent of re- 
cent cool dwarf-star photometry by Spitzer from Patten et al. ( 2006 ) 
and Kirkpat rick et al.| |20TTJ it is now clear that these models are 
inadequate. Specifically, while the largest observed colour over the 
spectral break of z — J ~ 3.5 (see Fig.[2| agrees with the AMES- 
dusty model, where z — J ~ 3, the 3.6/itn - 5.8/itn colour, which is 
proposed by Hsieh et al. (2012) as the best way to discriminate be- 
tween dwarf stars and galaxies, is underestimated by the model by 
over a magnitude. The results of Kirkpatrick et al7|( |201 lp in partic- 
ular, who observed mainly T-dwarf stars, show that dwarf stars can 
be significantly redder than predicted in the AMES-dusty model, 
with 3.6/im - 5.8/im ~ 0.5 perfectly possible. 

Despite these concerns, the 3.6^im - 5.8^im colour of TENIS- 
ZD1 reported by Hsieh et al. ( 2012) is still redder than observed in 
the coolest T-dwarf stars to date, albeit now by only ~ 2a. How- 
ever, it also exceeds the predictions of high-redshift galaxy tem- 
plates with reddening up to Ay ~ 1.0. Unfortunately, because 
|Hsieh et al. (2012) have neglected to provide the position of their 
putative high-redshift galaxy, we cannot check the IRAC photom- 
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Stacked photometry 
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Figure 8. The best SED fit (shown with a blue line) to the average photom- 
etry produced from a stack of our four most robust z ~ 7 galaxy candidates 
as shown in Fig.|4]and described in Section pTTj The best-fitting photomet- 
ric redshift for the stack is z = 6.98 ± 0.05, rising to z = 7.12 with 
allowance for a contribution from Lyman-a emission (see Table 2). The 
best fitting low-redshift model (z ~ 1.5) is shown as the red curve, which 
clearly cannot reproduce the large spectral break and flat near-IR colours of 
the stack photometry. 



etry, or indeed the photometry and claimed non-detections at any 
other wavelength, despite the fact that much of the key datasets 
are public. Our own SED fit to the photometry provided by |Hsieh| 
|et al.| J2012) supports the high-redshift solution (which they de- 
rived using the EAZY code; Brammer et al. 2008), but we note that 
the extreme redshift which results is largely a consequence of the 
arguably surprisingly bright magnitude reported at 5.8 fim, while 
exclusion of the low-redshift alternative solution depends critically 
on the error ascribed to the single J-band data point (in the ab- 
sence of any Y-band imaging of appropriate depth). Given that we 
cannot make our own photometric measurements for this object, it 
is difficult to comment further on the validity or otherwise of this 
proposed z = 7 — 8 galaxy. What we can say is that, given our 
own results, it would not seem unreasonable to detect one galaxy at 
z ~ 7 at J ~ 25 in the ECDFS, but also that it is somewhat unex- 
pected to see it detected in the two longest-wavelength IRAC bands, 
given that none of our own Ultra VISTA candidate objects are de- 
tected in the S-COSMOS 5.8 /im or 8 /J,m imaging. However, we 
note that the reported 5.8 /im and 8.0 pm magnitudes for TENIS- 
ZD1, of mAB ~ 23.3 and 23.6 respectively, are deeper than the 
S-COSMOS limits. 



the four individual galaxies averages to produce a relatively flat UV 
slope, as discussed further below. The derived photometric redshift 
for the stack is z = 6.98 ± 0.05, and the colour across the Lyman- 
break is z 1 — Y = 2.5 ± 0.2, consistent with the marginal detec- 
tion in the z'-band postage-stamp seen in Fig. [5] and yielding an 
extremely secure Lyman break. As detailed in Table 2, allowing 
the SED fitting to also include a contribution from Lyman-a emis- 
sion causes the estimated redshift of the stack to rise to z = 7.12, 
with the best fit implying a rest-frame Lyman-a equivalent width 
of£W ~40A. 



6.1.2 Extinction and star-formation rates 

The best-fitting SED to the four-object stack is based on a tau- 
model with a characteristic star-formation timescale r = 50 Myr, 
with solar metallicity and moderate reddening of Ay = 0.3 mag. 
This is somewhat lower than the value of Ay = 0.75 derived at 
z ~ 6 by Willott et alT] ( |2012) from a stack of 23 z ~ 6 galaxies 
with z' ~ 25.0. 

However, we caution that dust reddening is, of course, degen- 
erate with age/metallicity, and may also be exaggerated in single- 
component fits such as those utilised here, as there is evidence that 
a two-component star-formation history may better reproduce the 
combination of fairly blue UV continuum (as displayed by this 
stack - see below) and (relatively) red UV-to-optical colour seen 
in the SEDs of many z ~ 7 galaxies (Curtis-Lake et al., in prepa- 
ration). The uncertainty in derived Ay is important to bear in mind 
when considering the SFRs and sSFRs of these galaxies because, 
as can be seen from columns 12 and 13 in Table|3] the galaxies with 
the most extreme inferred SFRs all have correspondingly high val- 
ues of Ay. It is for this reason that we have also judged it helpful, 
in the final column of Table 2, to provide an alternative estimate of 
SFR, based simply on the Madau et al. ( 1998) conversion from UV 
flux density, assuming zero dust extinction. 

The SFR of the stack derived from the best-fitting SED as- 
suming a constant SFH is 98MQyr _1 or, if no dust extinction is 
assumed (i.e. using the Madau et al. 1998 prescription), the SFR = 
32 M yr _1 . From Table [3] it can be seen that, with zero dust, the 
SFRs of the individual ga axies lie in the range 25 — 50 Moyr -1 , 
but that adopting the best-fitting SED models can yield values as 
high as several hundred M0yr _1 , as a result of best-fitting values 
of Ay reaching as high as unity (although the single object with 
Ay > 1 is not from our 'robust' sample of objects). 



6 DISCUSSION 

6.1 Stacked photometry and physical properties 

6.1.1 Creation and analysis of the four-object stack 

An average stack was made from the postage-stamps of our four 
most robust z > 6.5 galaxy candidates as illustrated in Fig. [4] This 
was done both to create higher signal-to-noise photometry to better 
explore the implications of SED fitting, and also to double-check 
that no low-level i-band flux was present that might indicate con- 
tamination of our "robust" sample with objects at z < 6.5. Photo- 
metric measurements were made in the same way as for the indi- 
vidual candidates, using SEXTRACTOR in dual-image mode with 
selection in the Y+ J image. The resulting average photometry is 
included in Table [2] and the best-fitting SED is shown in Fig. [8] 
where the variation in the rest-frame UV photometry observed in 



6.1.3 UV slopes 

An estimate of the UV continuum spectral slope, /3 (fx oc A* 3 ), can 
be derived from the stacked near-infrared photometry in a number 
of ways. Following Bouwens et al. ( 2010 1 and Dunlop et al.| ( [2012] >, 
a simple estimate of ft can be derived from the J — H colour (delib- 
erately avoiding the Y-band which could be contaminated by Lya 
emission). This yields a highly uncertain value of /3 = — 1.1 ± 0.9, 
but the availability of the K s photometry enables a more accurate 
estimate based on power-law fitting to the J, H, K s photometry 
which gives /3 = —2.0 ± 0.2 (where the error is obtained from 
Monte Carlo simulations where the photometry is varied randomly 
according to the derived errors). Including the Y-band photome- 
try as well (i.e. fitting to Y, J, H, K a ) yields a slightly bluer value 
j3 = —2.16 ± 0.18, but this could potentially be affected by Lya 
emission and is, in any case, still clearly consistent with f3 ~ —2. 
We thus adopt /3 = — 2.0 ± 0.2 as the best estimate of /3 for the 
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four-object stack, and note that this is perfectly consistent with the 
value deduced by |Dunlop et al.H2012) and |McLure et al.|pOTT) 
for the most luminous of the z — 7 galaxies uncovered by the 
HST WFC3/IR surveys to date (albeit these are > 1.5 mag less 
luminous than those considered here). Our result is also consis- 
tent with the value obtained by Finkel stein et al.| |201 1 ) who found 
13 — — 2.04^Q 27 for a sample of lower-luminosity galaxies with 
L > 0.75L* at z = 7, with a weak trend to redder colours with in- 
creasing luminosity. However, our result is significantly lower than 
the value of /3 derived at z ~ 6 by Willott et al. (20121, who in- 
ferred a significantly redder value of /3 = —1.44 ± 0.10 from their 
stacked Y, J, H, K s photometry. 

6.1.4 Stellar masses and specific star-formation rates 

The stellar mass derived from the SED fitting of the stacked pho- 
tometry is M* = 4 x 10 9 M(-7), wh ich can be compared with the 



Willott et al. 



1 2012 1 of M* ~ 1O 1U M for bright 



value found by 

galaxies at z = 6 (although the results are consistent within the 
errors); we note that the galaxy detected by Hsieh et al. I2 0l2[ l, if 
really at z p hot ~ 7.8, has a substantially larger inferred stellar mass 
of M. = 3.2 x 1O 1O M . 

Based on the best-fitting model, the specific SFR (sSFR) of 
the stack is then ~ 30Gyr _1 , falling to ~ 8Gyr _1 if zero dust 
extinction is assumed. Clearly both these values (and indeed the 
values inferred from Table 2 for most of the individual objects) 
are higher than the average value of ~ 2 — 4Gyr _1 generally re- 
ported for studies of fainter Lyman-break galaxies at these redshifts 
(e.g. |Gonzalez et al. 20101, although they lie within the range of 
values found for individual objects by McLure et al. (201 1). How- 
ever it is probably premature to over-interpret these values as they 
could be biased high for potentially two reasons. First, if there re- 
ally is a significant range in sSFR at a given stellar mass, then since 
the UltraVISTA data has only just reached the depth required to se- 
lect genuine z ~ 7 galaxies, it is likely that those galaxies detected 
on the basis of their rest-frame UV flux will be biased towards high 
values of sSFR. Second, even if the intrinsic range in sSFR is small, 
photometric errors will conspire to yield a systematic over-estimate 
of the rest-frame UV luminosities of our candidates, a point which 
is also important for the estimation of number density as a function 
of magnitude (as discussed further below). 



6.2 Luminosity Function 

A full analysis of the implications of the new, luminous z ~ 7 
galaxies reported here for the form and evolution of the UV galaxy 
LF is deferred to a future paper (where we expect to be armed with a 
larger sample, and will have processed the required suite of detailed 
mock-data simulations). Nevertheless it is still instructive to briefly 
consider how our results compare with expectations based on ex- 
isting determinations of the galaxy UV LF at these early times. 
However, we caution that the LF parameters calculated below rely 
on an approximate estimate of the survey volume and should be 
considered a preliminary result. 

In the left-hand panel of Fig. 8 we illustrate how our results 
measure up against currently-published determinations of the UV 
galaxy LF at z ~ 7 as reported by |Ouchi et al. j2009), |McLure| 
[eTaLl ( |20l0) , |Oesch et al.|p0T0l > and |Bouwens et al.| ( |2011[ >. As 
can be seen from this plot, despite the continuing (important) de- 
bate over the faint-end slope (see Dunlop 2012), these indepen- 
dent determinations of the z ~ 7 LF are in generally good agree- 
ment, especially between the break luminosity at M*5oo — —20 



and Mi 500 ~ —21, the brightest luminosity effectively sampled by 
the pre-existing ground-based data (in fact the McLure et al. (2010) 
and Ouchi et al. (2009) LFs are so similar that only one can be use- 
fully plotted). This agreement also appears to extend out to the lu- 
minosity region probed here by UltraVISTA, at Mi 50 o < —22, but 
it is important to stress that this is largely a result of the common 
assumption of a Schechter form, and that these curves are simply 
extrapolations beyond the range of the pre-existing data. 

To locate our new results on this plot we have adopted the 
J-band magnitude of the stack (i.e. the average J-band magni- 
tude of our four most robust candidates), corrected to total magni- 
tude, to calculate the rest-frame UV absolute magnitude at 1500A, 
Misoo = —22.5. We have then conducted a set of simulations and 
find that our selection criteria, based on the Y+ J image, implies 
that incompleteness and flux-boosting virtually cancel each other 
in this magnitude regime, and that a relatively modest flux-boosting 
correction of only ~ 0.1 mag is merited. This leads to a data-point 
in Fig. 8 plotted at A/1500 = —22.4, at a comoving number den- 
sity based on the discovery of four galaxies within our estimated 
comoving survey volume of 5 x 10 6 Mpc 3 , based on our effective 
survey area of 0.9 deg 2 (excluding the area unavailable for galaxy 
selection due to diffraction halos/spikes from bright stars) and as- 
suming that our selection function is effectively unity over the red- 
shift range z ~ 6.5 — 7.2. The size of the bin in M1500 was chosen 
to be 1 mag to indicate the range of the observed J-band measure- 
ments for our four most secure galaxies (as seen in Table[2](, while 
the error bar in density is based purely on the Poisson noise. Finally, 
to indicate possible extremes, we also show where this data-point 
would lie in the unlikely circumstance that all ten of the objects 
presented here are truly at z > 6.5, or alternatively, if only one 
high-redshift galaxy were to be confirmed. 

It is clear from this plot that, while our results are much less 
extreme than would be implied from the results of Capak et al. 
(2011), nevertheless our discovery of four robust z ~ 7 galaxies 
at M1500 — —22.4 yields a somewhat higher number density than 
would be anticipated on the basis of simply extrapolating the ex- 
isting LFs. For example, the McLure et al. (2010) LF predicts only 
~ 0.8 galaxies with J < 25 in our survey volume. 

Given the small-number statistics, and the potential for cosmic 
variance (although with the low sample size available here, Poisson 
noise dominates that from cosmic variance), it would be premature 
to suggest that our results contradict the steep exponential fall-off at 
the bright end of the LF implied by the current Schechter-function 
fits to the z ~ 7 LF. However, it also important to stress that, as is 
clear from Fig. 8, our results are certainly not in conflict with the 
existing data at this redshift, and would in fact lie naturally on a 
power-law extrapolation from the data at M* 50 o — —21. To illus- 
trate this we have plotted a simple double power-law LF through 
the data, with the functional form: 



4>(M) 



1Q0.4(A+l)(Af-A/*) _|_ lQ0.4(B+l)(Af-Af*) ' 



(1) 



We set cf>* , M* and the faint end slope, A (= a), to the values 
determined in the Schechter-function fit by Mc Lure et al.| l [20To) l, 
and have simply set the value of the bright-end slope, B, to make 
the function pass through our data-point (B = —4.2). Even without 
a detailed analysis it is obvious that this functional form provides a 
good description of all the available data, indicating that it may be 
premature to conclude that the LF at 2 ~ 7 is best described by a 
Schechter function. 

Finally, in the second panel of Fig. 8, we show the same z ~ 7 
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Figure 9. The left-hand panel shows the 2 = 7 galaxy UV (~ 1500 A) LF with data shown from McLure et al. (2010 1, Bouwens et al 



2010a|b) and |Ouchi et al.|{2009) overplotted on the best-fitting Schechter functions derived by Bouwens et al. 



j201l) , |Castellano etal. 



2011) , |Oeschetal.| l 20 1 1 and McLure et al. 



2010) (note the LFs derived by Ouchi et al. (2009) and McLure et al. (2010) are essentially identical). The estimated data-point from our new UltraVISTA 



study is shown in red, with Poissonian errors on the number density, with the upper, central and lower point given by ten, four or only one candidate in our 
sample being a confirmed as a z > 6.5 galaxy. A double power-law curve is also included for comparison, where the faint-end slope is matched to that 
of McLure et al. (2011) and at the bright end we adopt B= —4.2 (see Section[6]for p arameterisation). In the right-hand plot we show the same data, but this 
time compare with the the z = 6, 7, 8 LFs from McLure et al. 1 2009 1, McLure et al. ( 2010 1 and Bradley et al. 1 2012 1 respectively (see text for further details 
and discussion). 



data points, but here we plot the best-fitting galaxy UV LFs at 
z ~ 6 (McLure et al. 2009; Willott et al 2012), z ~ 7 (McLure 
et al. 2010) and z ~ 8 (Bradley et al. 2012). This figure demon- 
strates that there is in fact currently very little (if any) significant 
evidence for evolution in the bright end of the galaxy LF in the 
redshift range z ~ 6 — 8. We note that the only other work directly 
probing the galax y LF at Mi 500 < —22 at these epochs is the 
study at z ~ 6 by |Willott et al. 1 2012 1, who exploited ~ 4 deg 2 of 
ground-based CFHT imaging to search for 2 ~ 6 galaxies brighter 
thanMi35o ~ -21.5. |Willott et al.| l [20T2] > found 40 galaxies, many 
with M1350 < -22, but relatively few with Mi 350 < -22.5. 
Thus, even at z ~ 6 the basic statistics at Muv < —22.5 re- 
main poor, but nevertheless a maximum likelihood fit to the 2 ~ 6 
LF led Wi llott et al.| j2012) to conclude in favour of the steep ex- 
ponential cutoff implied by the Schechter-function fit of McLure 
et al. (2009). However, the near-infrared data analysed by Willott 
et al. (2012) are somewhat inhomogeneous, of inadequate depth to 
detect the 2 ~ 6 galaxies in J-band in two of the four survey fields, 
and lack the F-band data which, ideally, would enable determina- 
tion of M1500 free from corrections due to IGM extinction (which 
can impact on the 2'-band at these redshifts). Moreover, in the sam- 
ple reported by Willott et al. (2012) the COSMOS field alone con- 
tains two 2 ~ 6 galaxies with J-band magnitudes brighter than 
the most luminous galaxy candidate (as based on the 2-band mag) 
reported from the full 4 deg 2 survey. We conclude that the form of 
the bright-end of the LF is still open to debate at both 2 ~ 6 and 
2 ~ 7. At 2 = 8, Bra dley et al-H2012) have also claimed to find no 
evidence for an excess of sources at the bright-end. However, due 
to the (relatively) small effective field size of the current BoRG sur- 



vey, the galaxies detected are at L < 2L* (the galaxies presented 
here are ~ 9L*) and so the constraints on the very bright-end of the 
2 ~ 8 LF remain weak. The final UltraVISTA imaging is designed 
to improve our knowledge of the bright end of the galaxy UV LF 
at z ~ 8 as well as at 2 ~ 7. 



7 CONCLUSION 

We have exploited the new, deep, near-infrared UltraVISTA imag- 
ing of the COSMOS field, in tandem with deep optical and mid- 
infrared imaging, to conduct a new search for luminous galaxies 
at redshifts 2 ~ 7. The unique multi-wavelength dataset provided 
by VISTA, CFHT, Subaru, HST and Spitzer over a common area 
of ~ 1 deg 2 has allowed us to select galaxy candidates at red- 
shifts 2 > 6.5 by searching first for UltraVISTA Y+ J-detected 
(< 25 mag) objects which are undetected in the CFHT and HST 
optical data. 

This sample was then refined using a photometric redshift fit- 
ting code, enabling the rejection of lower-redshift galaxy contami- 
nants, and cool galactic M, L, T dwarf stars; brown-dwarf contami- 
nation is a much more serious problem for wide-area ground-based 
surveys than for deeper/narrower HST WFC3/IR surveys, and so 
we have taken great care to utilise the full multi-wavelength dataset 
(fitting the latest stellar templates), including IRAC colours, to min- 
imise dwarf-star contamination. 

The final result of this process is a small sample of (at most) 10 
credible galaxy candidates at 2 > 6.5 (from over 200,000 galaxies 
detected in the year-one UltraVISTA data). The first four of these 
objects appear to be robust galaxies at 2 > 6.5, and fitting to their 
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stacked spectral energy distribution yields z p hot = 6.98 ± 0.05 
with a stellar mass M* ~ 5 x 10 9 Mq and rest-frame UV spectral 
slope /3 ~ —2.0 ± 0.2 (where fx oc A 8 ). The next three are also 
good candidates for z > 6.5 galaxies, but the possibility that they 
are low redshift interlopers or dwarf stars cannot be excluded. Our 
final subset of three additional candidates is afflicted not only by 
potential dwarf-star contamination, but also contains objects likely 
to lie at redshifts just below z = 6.5. 

We have also been able to demonstrate that the three even- 
brighter z > 7 galaxy candidates reported in the COSMOS field 
by |Capak et al.| ( |2011[ ) are in fact all lower-redshift galaxies at 
z ~ 1.5 — 3.5. Consequently the new z ~ 7 galaxies reported here 
are the first credible z ~ 7 Lyman-break galaxies discovered in the 
COSMOS field and, as the most UV-luminous discovered to date at 
these redshifts, are prime targets for deep follow-up spectroscopy. 
We have investigated the physical properties of these galaxies as 
inferred from the broad-band photometry, and have discussed the 
uncertainties in, and implications of their estimated star-formation 
rates and stellar masses. Finally, we have considered only briefly 
the implications of the inferred number density of these "bright" 
galaxies for the form of the galaxy luminosity function at these 
early epochs, deferring a full re-analysis of the z ~ 7 UV luminos- 
ity function to a future paper (when we expect to be armed with a 
larger sample of objects). 
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